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The excitation function for the pp — > ppK + K~ reaction revealed a significant enhancement close 
to threshold which may plausibly be assigned to the influence of the pK~ and K + K~ final state 
interactions. In an improved reanalysis of COSY-11 data for the pp — ► ppK + K~ reaction at ex- 
cess energies of Q = 10 MeV and 28 MeV including the proton- K~ interaction the enhancement is 
confirmed. Invariant mass distributions for the two- and three-particle subsystems allow to test at 
low excess energies the ansatz and parameters for the description of the interaction in the ppK + K~ 
system as derived from the COSY-ANKE data. Finally, based for the first time on the low energy 
K + K~ invariant mass distributions and the generalized Dalitz plot analysis, we estimate the scat- 
tering length for the K + K~ interaction to be \FLe(a K + K -)\ — 0.5^'° fm and Im(a K + K -) = 3.0 ± 
3.0 fm. 
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I. INTRODUCTION 

The basic motivation for investigating the pp — > 
ppK + K~ reaction near the kinematical threshold is com- 
prehensively reviewed in [l[ , as an attempt to understand 
the nature of the scalar resonances / (980) and a (980), 
whose masses are very close to the sum of the K + and 
K~ masses. Besides the standard interpretation as qq 
mesons 0, these resonances were also proposed to be 
qqqq states @, KK molecules 0, [f|, hybrid gq/meson- 
meson systems @ or even quark- less gluonic hadrons 0] • 
The strength of the KK interaction is a crucial quan- 
tity regarding the formation of a KK molecule, whereas 
the KN interaction is of importance in view of the vig- 
orous discussion concerning the structure of the excited 
hyperon A(1405) which is considered as a three quark 
system or as a KN molecular state @. Additionally, 
these interactions appear to be very important also with 
respect to other phenomena, like possible kaon conden- 
sation in neutron stars [9| , or production of strange par- 
ticles immersed in a dense nuclear medium studied by 
means of heavy ion collisions [IE (HI EH EH ■ 

Measurements of the near threshold pp — * ppK + K~ 
reaction have been made possible by beams of low emit- 
tance and small momentum spread available at storage 
ring facilities, in particular at the cooler synchrotron 
COSY at the research center in Julich in Germany [ijj]. 
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A precise determination of the collision energy, in the or- 
der of fractions of MeV, permitted to deal with the rapid 
growth of cross sections [l5| and thus to take advantage 
of the threshold kinematics like, e.g., full phase space 
coverage achievable with dipole magnetic spectrometers 
being rather limited in geometrical acceptance. Early 
experiments on K + K~ pair production at COSY con- 
ducted by the COSY-11 collaboration revealed, however, 
that the total cross section at threshold is by more than 
seven orders of magnitude smaller than the total proton- 
proton production cross section making the study diffi- 
cult due to low statistics [IE E3, EH ■ A possible influence 
from the /o or ao mesons on the K + K~ pair production 
appeared to be too weak to be distinguished from the di- 
rect production of these mesons based on the COSY-11 
data [ri| . Recent results obtained by the ANKE collab- 
oration with much higher statistics can also be explained 
without the need of including the scalars /o or ao |19ll20j. 
However, the combined systematic collection of data be- 
low HE SEES] and above EE El the (j) meson threshold 
reveal a significant signal in the shape of the excitation 
function which may be due to the K~p and perhaps also 
to the K + K~ interaction. This signal is based on the 
COSY-11 data which, as indicated by authors of arti- 
cle EH , were analyzed calculating the acceptance without 
the inclusion of the pK~ interaction. Therefore a more 
detailed analysis of the COSY-11 data at excess energies 
of Q = 10 MeV and 28 MeV including now studies of 
both the differential cross section distributions and the 
strength of the final state interaction between the K + 
and K~ mesons was performed. The analysis is based 
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on a generalization of the Dalitz plot for four particles as 
proposed by Goldhaber et al. [13, HH [H| • 



II. EXCITATION FUNCTION FOR THE NEAR 
THRESHOLD pp -> ppK+R- REACTION 

The measurements of the pp — > ppK + K~ reaction 
were conducted at low excess e nerg ies by the collabora- 
tions ANKE [H|, COSY-11 HE Eali3 and DISTO 
The achieved results are presented in Fig. [T] together 
with curves representing three different theoretical ex- 
pectations normalized to the DISTO data point at 
Q = 114 MeV. The dashed curve represents the energy 
dependence from four-body phase space when we assume 
that there is no interaction between particles in the final 
state. These calculations differ from the experimental 
data by two orders of magnitude at Q = 10 MeV and 
by a factor of about five at Q = 28 MeV. Hence, it is 
obvious, that effects of final state interactions cannot be 
neglected in the ppK + K~ system (25[. Inclusion of the 
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FIG. 1: (Color online) Excitation function for the pp — > 
ppK + K~ reaction. Triangle and circles represent the DISTO 
and ANKE measurements, respectively. The four points close 
to the threshold are results from the COSY-11 measurements. 
The curves are described in the text. 

pp-FSl (dashed-dotted line in Fig. [T]) , by folding its pa- 
rameterization known from the three body final state [26[ 
with the four body phase space, is closer to the exper- 
imental data, but does not fully account for the differ- 
ence [lH . The enhancement may be due to the influence 
of pK and K + K~ interaction which was neglected in 
the calculations. Indeed, as shown by authors of refer- 
ence [l^ll^l the inclusion of the pK~-FSl (solid line) re- 
produces the experimental data for excess energies down 



to Q = 28 MeV. These calculations of the cross section 
were accomplished under the assumption that the over- 
all enhancement factor, originating from final state in- 
teraction in the ppK + K~ system, can be factorised into 
enhancements in the pp and two pK~ subsystems (l9| : 

Ffsi = F pp (q) ■ F piK - (fe) • F P2K - (fc 2 ) , (1) 

where fci, ki and q stand for the relative momenta of 
the particles in the first pK~ subsystem, second pK~ 
subsystem and pp subsystem, respectively. The factors 
describing the enhancement originating from the pK~ 
FSI are parametrized using the scattering length approx- 
imation. It is important to note that the inclusion of 
the pp and pK~ final state interaction is not sufficient 
to describe the data very close to threshold (see Fig.[T]). 
This enhancement may be due to the influence of the 
K + K~ interaction, which was neglected in the calcula- 
tions jyj. However, as pointed out in Ref. the ob- 
served increase of the total cross section near threshold 
may be due to the neglect of the pK~-FSl in the cal- 
culations of the COSY-11 acceptance. As a consequence 
the extracted cross sections would decrease, if this inter- 
action was taken into account during the analysis of the 
experimental data. This concern encouraged us to check 
quantitatively the influence of the interaction in the pK~ 
subsystem on the acceptance of the detection setup. In 
addition, absolute values for the differential distributions 
of the pK and ppK invariant masses were extracted and 
generalized Dalitz plot analysis of the data in view of the 
K + K~ interaction, was performed. 



III. MEASUREMENTS OF THE pp -» ppK+K' 
REACTION PERFORMED WITH THE COSY-11 
MAGNETIC SPECTROMETER 

The measurements of the pp — > ppK + K~ reaction 
close to threshold have been conducted using the cooler 
synchrotron COSY and the COSY-11 detector sys- 
tem [13] shown schematically in Fig.O The target, being 
a beam of H 2 molecules grouped to clusters of up to 10 5 
atoms .28] , crosses perpendicularly the proton beam cir- 
culating in the ring. If a collision of protons leads to the 
production of a K + K~ meson pair, then the reaction 
products, having smaller momenta than the circulating 
beam, are directed by the magnetic dipole field towards 
the COSY-11 detection system and leave the vacuum 
chamber through a thin exit foil [2?| . Tracks of positively 
charged particles, registered by drift chambers, are traced 
back through the magnetic field to the nominal interac- 
tion point leading to a momentum determination. The 
knowledge of the momentum combined with an indepen- 
dent measurement of the velocity, performed by means 
of the scintillation detectors SI and S3, permits to iden- 
tify the registered particles and to determine their four 
momentum vectors. Knowing both the four momenta 
of the three positively charged ejectiles and the proton 
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FIG. 2: Schematic view of the COSY-11 detector with a typ- 
ical event of the pp — > ppK + K~ reaction channel. For the 
description see text. 



beam momentum one can calculate the mass of the un- 
observed system X~. In the case of the pp — > ppK + K~ 
reaction this should correspond to the mass of the K~ 
meson, but we observe also background originating partly 
from the pp — > ppn + X~ reaction, where the n + was 
misidentified as a K + meson, and in part due to the 
K + meson being associated with the hyperons A(1405) 
or £(1385) production [13, EE]- This background, how- 
ever, can be completely removed by demanding a signal 
in the silicon pad detectors (mounted inside the dipolc) 
at the position where the K~ meson originating from 
the pp — > ppK + K~ reaction is expected. This clear 
identification allows to select the pp — > ppK + K~ events 
and to determine the total and differential cross sections. 
A more detailed description of the experiment and data 
evaluation is given in LLZI, [lH, [27|, [29| . 



IV. DIFFERENTIAL OBSERVABLES FOR 
COSY-11 DATA AT Q = 10 MeV AND Q = 28 
MeV 

In order to check the sensitivity of the result to the 
assumption of the pK~ final state interaction we derived 
the distributions of the differential cross section assum- 
ing that the acceptance depends only on the pp-FSI [4|| . 
Then we calculated the acceptance with inclusion of the 
pp- and pK~-FSI, and derived analogous distributions. 
In this calculations we assumed the factorisation of the 
final state interaction given by Eq. Q] and used the pK~ 
scattering length a p K- = (0 + 1.5i) fin The results 
are presented in Fig. [3] for data at Q = 10 MeV and in 
Fig. 0] for Q = 28 MeV. The distributions obtained under 



TABLE I: The differential cross sections of the pp — > 
ppK + K~ reaction as a function of invariant masses for differ- 
ent subsystems. The values of My denote the center of the 
invariant mass bins of 2.5 MeV/c 2 and 7.0 MeV/c 2 width for 
Q = 10 MeV and Q = 28 MeV, respectively. 



Q = 10 MeV 



M PP (^) 



1.8778 
1.8803 
1.8828 
1.8853 



ipp \ GcV/c' A 



216 ± 53 
106 ± 34 
43 ± 22 
14 ± 14 



M, 



pK+ 



K + ^GcV/c 



1.4332 
1.4357 
1.4382 
1.4407 



107 ± 34 
122 ± 39 
125 ± 42 
32 ± 21 



AI 



pK- 



■ GeV ■ 



_ V GeV/c 2 



1.4332 
1.4357 
1.4382 
1.4407 



173 ± 47 
145 ± 42 
36 ± 21 
13 ± 11 



M 



K+K- 



M 



0.9887 
0.9912 
0.9937 
0.9962 

UcV \ 



dM K+K- VGcV / 



169 ± 44 
174 ± 51 
35 ± 21 
± 9 



K+ V G eV/c^ 



AI, 



2.3715 
2.3740 
2.3765 
2.3790 

ppK~ 



± 13 
68 ± 28 
164 ± 43 
99 ± 38 



(^) 



~n5~ 



GeV/c 



2.3715 
2.3740 
2.3765 
2.3790 



60 ± 30 
115 ± 39 
127 ± 39 
87 ± 31 



Q = 28 MeV 



Mpp (W) 



1.880 
1.887 
1.894 
1.901 



388 
346 
148 
48 



> GeV/c 2 

± 92 
± 87 
± 67 
± 48 

I nb 



1.435 
1.442 
1.449 
1.456 



255 
261 
287 
63 



AI, 



pK- 



GeV/c 

± 72 
± 76 
± 83 
± 37 

I nb 



AI 



K+K 



1.435 
1.442 
1.449 
1.456 



581 
135 
97 
25 



^GeV/c 2 
± 117 

± 56 
± 46 
± 21 

nb 



0.991 
0.998 
1.005 
1.012 



M K+K - 



221 
454 
230 
38 



AI 



pK + 



AI, 



ppK 



2.374 
2.381 
2.388 
2.395 



20 
61 
247 
566 



2.374 
2.381 
2.388 
2.395 



204 
284 
176 
216 



( GcV/c 2 > 

± 70 
± 114 
± 70 
± 22 

{ nb \ 

\ GcV/c 2 I 

± 20 
± 36 
± 72 
± 121 

(g ' 
± 68 
± 79 
± 63 
± 69 



eV/c 



both assumptions are almost identical, which shows that 
the acceptance of the COSY-11 detection setup is only 
very weakly sensitive to the interaction between K~ and 
protons. Thus, the observed enhancement in the excita- 
tion function cannot be explained by approximations in 
the determination of the detection efficiency as suspected 
by p^ |. This justifies the assumption made in the orig- 
inal analysis, where the efficiency was calculated taking 
into account the pp-FSl only. 

The derived values of differential cross sections are listed 
in Table [U This result constitutes an additional informa- 
tion to the total cross sections published previously [lal . 
The values of the cross sections in the former analysis [18[ 
were determined using the total number of events iden- 
tified as a pp — > ppK + K~ reaction and the total ac- 
ceptance of the COSY-11. Now after the determination 
of the absolute values for the differential distributions 
one can calculate the total cross sections in a less model 
dependent manner regardless of the assumption of the 
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FIG. 3: (Color online) Differential cross sections for the 
pp — > ppK + K~ reaction at Q = 10 MeV. Circles and dashed 
bars (in black) denote spectra where the acceptance was deter- 
mined taking into account only the pp-FSI, and triangles with 
solid bars (in red) denote results where pp- and pK~-FSl was 
taken into account in the acceptance calculations. They are 
hardly distinguishable. Vertical bars indicate statistical er- 
rors, whereas the horizontal bars stand for the invariant mass 
intervals for which the cross section values were established. 



pp-FSl. The cross sections, calculated for both excess 
energies as an integral of the M pp distribution derived 



with the inclusion of the pK 
culations: 



-FSI in the acceptance cal- 



Otot 



da 



amount to Ctot = (0.95 ± 0.17) nb for measurement at 
Q = 10 MeV and cr tot = (6.5 ± 1.1) nb for Q = 28 MeV. 
These results are statistically consistent with the previ- 
ously evaluated total cross sections. However, the values 
extracted in the actual analysis with smaller error bars 
are larger by about 20% for Q = 10 MeV and 50% for 
Q = 28 MeV, which strengthens the confidence to the 
observed enhancement at threshold. 
The determination of the absolute values for the differen- 
tial cross sections permitted us to establish the absolute 



FIG. 4: (Color online) Differential cross sections for the pp — * 
ppK + K~ reaction at Q = 28 MeV. For the description see 
caption of Fig. [3] 



values for the following ratios at the close to threshold 



region [46 



RpK 
RppK 



da/dM pK - 
da/dM p K+ 
da/dM ppK ~ 
da/dM ppK + 



If pK + and pK~ interactions were the same, the distri- 
bution of R p k as well as R pp k should be flat and equal 
to unity. But as one can see in Fig. [5] and as presented 
already in the previous publication by COSY-11 [l8[ and 
ANKE [19] for both excess energies R p k is far from being 
constant and increases towards the lower M p x invariant 
masses. This effect might be due to the influence of the 
pK~ final state interaction. Similarly the distributions 
of R pp k differs from expectations assuming only FSI in 
the pp system. This is a confirmation of effects observed 
also by the ANKE collaboration at higher excess ener- 
gies [131 • The determination of the absolute values of the 
ratios R p k and R pp k allows to compare the parameters 
of the scattering length a pK - derived by the ANKE group 
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FIG. 5: (Color online) The distributions of the absolute values 
of ratios R p k and R vv k for data at Q = 10 MeV ((a) and 
(b)) and Q = 28 MeV ((c) and (d)). Solid curves represent 
theoretical expectations calculated taking into account pp and 
pK~ final state interaction. It should be stressed that the 
absolute values are not fitted, but they are result of the model 
and used parameters 



Generalization of the Dalitz plot: 
approach 



Goldhaber 



There are many different types of generalizations of 
the Dalitz plot for four-body final states [H], [3f| . Here 
a g eneralization proposed by Goldhaber is presented [22], 

MM- 

Consider a reaction like a + b — > 1 + 2 + 3 + 4 yield- 
ing four particles with masses m, and total energy yfs 
in the centre-of-mass frame. The probability, that the 
momentum of the i th particle is in a range d 3 pi is given 
by: 



1 



d LZ P = d i p 1 d A p 2 d i p ? ,d i p 4 - 

VI P. Pi P4 lQ ElE2EsE4 

><^(epjJmE^-^ji m i 2 - ( 2 ) 



where Ei = \/pf + mf denotes the energy of the i th par- 
ticle (c = 1) and M denotes the invariant matrix element 
for the process. Assuming that the matrix element M de- 
pends only on the invariant masses of the two- and three 
particle subsystems [34| , the distribution given by Eq. [2] 
can be expressed in some choice of five independent in- 
variant masses (47| . An especially convenient choice is 
Ml 2, M^Ml^Ml 24, and Mf 34 Using such vari- 

ables, after integrations, one obtains an event distribu- 
tion in the following form: 



(from data at excess energies above the (f> meson pro- 
duction threshold), to the present data near the K + K~ 
threshold. As demonstrated in Fig. [5J simulations tak- 
ing into account the pK~ final state interaction with the 
scattering length a pK - = (0 + 1.5i) fm determined by 
the ANKE group from data at significantly higher excess 
energies reproduce very well the distributions of R p k and 
RppK near the threshold. 



V. ANALYSIS OF THE K + R- FINAL STATE 
INTERACTION 

A factorization ansatz for the pp and pK~ interac- 
tion underestimate the excitation function for the pp — > 
ppK + K~ reaction very close to threshold indicating that 
in this energy region the influence of the K + K~ interac- 
tion is significant and cannot be neglected. The interac- 
tion may manifest itself even stronger in the distributions 
of the differential cross sections jl5|. This observation 
demands an analysis of the double differential cross sec- 
tions for the low energy data at Q = 10 MeV (27 events) 
and Q = 28 MeV (30 events), in spite of the quite low 
statistics available. 



d 5 P = — \M\ 2 —= dMi 2 dMl 4 dMl 4 dMl 24 dMl 34 , 
os y — 1> 

(3) 

where B is a function of the invariant masses with the 
exact form to be found in Nyborg's work (34|. An advan- 
tage of the choice of mass variables used here is the high 
symmetry of the function B (34[, which is very usefull 
in the consideration of the boundary of the kinemati- 
cally allowed region in the (Mf 2 , M| 4 , Mf 4 , Mf 24 , Mf 34 ) 
space defined by B = 0. If we additionally change the 
integration variables to invariant masses the projection 
of the physical region on the (M12, M3 4 )-plane gives a 
convenient and simple distribution, which can be used 
to analyse the final state interaction in the same way as 
in case of three particles. Such an anal ysis was origi- 
nally made by Goldhaber et al. in 1963 [22j, |23j and is 
called Goldhaber plot. As it is shown in Fig. [S] the kine- 
matically allowed region in the Goldhaber plot is a right 
isosceles triangle within which the area (contrary to the 
Dalitz plot) is not proportional to the phase space vol- 
ume [34|. Moreover, even in the case of the absence of 
any final state interaction the density of events on the 
Goldhaber plot is not homogeneous. However, it is still 
a very convenient tool due to its lorentz invariance and 
simple boundary equations [34| : 

M12+M34 = Vs , M\i — mi+m 2 , and M 34 = m 3 + m 4 . 
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FIG. 6: Goldhaber plots for the pp — ► ppK + K~ reaction. 
The solid lines of the triangles show the kinematically allowed 
boundaries. Raw data are shown in Figs, (a) and (b) as 
black points. The superimposed squares represent the same 
distributions but binned into intervals of AM = 2.5 MeV/c 2 
(AM = 7 MeV/c 2 ) widths for an excess energy of Q = 10 MeV 
(28 MeV), respectively. The size of the square is proportional 
to the number of entries in a given interval. In Figs, (c) 
and (d) Monte Carlo results are presented. In the simulated 
distributions both the pp and the pK~-FSl are taken into 
account. 



B. Determination of the K + K scattering length 

Complementary to previous derivations [H, H3, [H, [39[ 
here we estimate the K + K~ scattering length directly 
from the low energy differential mass distributions of 
K + K~ and pp pairs from the ppK + K~ system pro- 
duced at threshold. The raw data (represented by black 
points in Figs. Ufa) and M,b)) were first binned into in- 
tervals of AM = 2.5 MeV/c 2 width for the measurement 
at Q = 10 MeV and intervals of AM = 7 MeV/c 2 for the 
data at Q — 28 MeV, and then for each bin corrected for 
the acceptance and detection efficiency of the COSY-11 
facility [4(| . The resulting Goldhaber plots are presented 
together with the raw distributions in Figs. [Ha) and[5fb). 
Figures [6fc) and Efd) show corresponding distributions 
simulated with Monte Carlo method taking into account 
the pp and pK~ interaction according to the factorization 
ansatz [19 ] . 

In order to estimate the strength of the K + K~ interac- 
tion, the derived cross sections were compared to results 
of simulations generated with various parameters of the 
K + K~ interaction taking into account strong final state 
interaction in the pp and pK~ subsystems. To describe 
the experimental data in terms of final state interactions 
between i) the two protons, ii) the K~ and protons and 
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FIG. 7: x 2 ~ Xmin distribution as a function of \Re{a K + K - 



Xmin 



denotes the absolute minimum with 



and Im(a K + 1 

respect to parameters a, lit^a^^- )|, and Im(a K + K -). 

iii) the K + and K~ , the K + K~ enhancement factor was 
introduced such that Eq. [T] changes to: 



*FSI 



F PP (q)-F piK - {k x )-F V2K - (k 2 )-F K+K - (k 3 ) . (4) 



As for the case of the pK -FSI, the F K + K - was calcu- 
lated in the scattering length approximation: 



F, 



1 



K+K- — 



1 - i k 3 a K + K - 



(•5) 



where cik+k- is the effective K + K~ scattering length 
and &3 stands for the relative momentum of the kaons 
in their rest frame. Using this parametrization we com- 
pared the experimental event distributions to the results 
of Monte Carlo simulations treating the K + K~ scatter- 
ing length as an unknown parameter, which has to be 
determined. In order to estimate the real and imagi- 
nary part of a K+K- we constructed the Poisson likeli- 
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hood x 2 statistic derived from the maximum likelihood 
method 

x ^a K+K - 1 a) = 2-Y J WNt-Nt + N^H-^)] , (6) 

i ' 1 

where Nf denotes the number of events in the i th bin 
of the experimental Goldhaber plot, Nf stands for the 
content of the same bin in the simulated distributions, 
and a is the normalization factor. The data collected at 
both excess energies have been analysed simultaneously. 
The obtained x 2 distributions (suppressed by its mini- 
mum value) as a function of the real and imaginary part 
of the K + K~ scattering length are presented in Fig. [7] 
The best fit to the experimental data corresponds to 
\Re(ax+K-)\ = 0.5 *q 5 fm and Im(aK+K- ) = 3 ± 3 fm. 
The final state interaction enhancement factor F K + K - in 
the scattering length approximation is symmetrical with 
respect to the sign of Re(ax+K-)i therefore only its ab- 
solute value can be determined. 



VI. CONCLUSIONS 

In conclusion, the more detailed analysis of the COSY- 
11 data with inclusion of the pK~ final state interaction 
did not change significantly the result of the previous 
analysis Moreover, the new more precise determi- 

nation of the total cross sections from the differential M pp 
distributions even increased the enhancement at thresh- 
old. 

In addition the analysis of the pp — > ppK + K~ reaction 



has been extended to the determination of differential 
cross sections in view of the K + K~ final state interac- 
tion. The extracted K + K~ scattering length amounts 
to: 

a K +K- = [(0.5 tU) + i(3 ± 3)] fm. 

Due to the low statistics the uncertainties are rather 
large. In this analysis we cannot distinguish between 
the isospin 1 = and 1=1 states of the K + K~ system. 
However, as pointed out in [43| . the production with 1 = 
is dominant in the pp — * ppK + K~ reaction independent 
of the exact values of the scattering lengths. 

Regarding the comparison of the interactions in the 
pK~ , pK + , ppK~ and ppK + subsystems, the absolute 
ratios determined from the COSY-11 data measured at 
Q = 10 MeV and Q = 28 MeV are consistent with the 
predictions based on the parametrization introduced in 
reference [l9| and on the values of the scattering length 
a p K- extracted from the ANKE data at higher excess 
energies fl9j . 
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